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Abstract 

Equations that follow from the Navier-Stokes equation and incom- 
pressibility but with no other approximations are "exact.". Exact equa- 
tions relating second- and third-order structure functions are studied, as 
is an exact incompressibility condition on the second-order velocity struc- 
ture function. Opportunities for investigations using these equations are 
discussed. Precisely defined averaging operations are required to obtain 
exact averaged equations. Ensemble, temporal, and spatial averages are 
all considered because they produce different statistical equations and 
because they apply to theoretical purposes, experiment, and numerical 
simulation of turbulence. Particularly simple exact equations are ob- 
tained for the following cases: i) the trace of the structure functions, ii) 
DNS that has periodic boundary conditions, and iii) an average over a 
sphere in r-space. The last case (iii) introduces the average over orienta- 
tions of r into the structure function equations. The energy dissipation 
rate e appears in the exact trace equation without averaging, whereas in 
previous formulations e appears after averaging and use of local isotropy. 
The trace mitigates the effect of anisotropy in the equations, thereby re- 
vealing that the trace of the third-order structure function is expected 
to be superior for quantifying asymptotic scaling laws. The orientation 
average has the same property. 

1 Introduction 

Equations relating statistics for turbulence studies, such as Kolmogorov's 
(1941) equation are asymptotic equations. This has required experimenters 
to seek turbulence that satisfies the criteria of the asymptotic state. The 
present approach is to derive exact statistical equations. These can be used to 
determine all effects contributing to the balance of statistical equations. By 
"exact" we mean that the equations follow from the Navier-Stokes equation and 
the incompressibility condition with no additional approximations. 

Exact equations have the potential to detect the limitations of direct 
numerical simulation (DNS) and of experiments and to study the approach to 
local homogeneity and local isotropy and scaling laws (Hill 2001). For those 
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purposes, the averaging operation must be exactly denned and implemented; 
that is done here. The methods developed here can be used on the exact 
structure-function equations of all orders N; those equations are in Hill (2001). 
It is useful to further investigate the exact second-order (N — 2) equation, 
which relates second- and third-order structure functions, because it has special 
simplifications that the higher-order equations (N > 2) do not possess and 
because the second-order equation is the most familiar. Exact equations satisfy 
the need perceived by Yaglom (1998) for careful derivation of dynamic-theory 
equations and the perceived value placed by Sreenivasan & Antonia (1997) on 
aspects of turbulence that can be understood precisely. Experimental data have 
been used to evaluate the balance of Kolmogorov's equation (Antonia, Chambers 
& Browne 1983; Chambers & Antonia 1984) and generalizations of it (Lindborg 
1999; Danaila et al. 1999 a,b; Antonia et al. 2000). This report supports such 
experimental work as well as precise use of DNS by giving exact equations to 
be used in such evaluations. The connection between the derivations presented 
here and any experiment or DNS is important because the equations relate 
several statistics and therefore are most revealing when data are substituted 
into them. 

The equations derived here are exact for every flow, whether laminar or 
turbulent, provided that no forces act on the fluid at the points of measurement, 
which points are denoted below by x and x'. Forces can be applied near the 
point of measurement; e.g., the equations are exact for hot-wire anemometer 
supports just downstream of the measurement points. The cases of forces 
at the points of measurement and throughout the fluid are considered in Hill 
(2002), wherein the Kolmogorov flow computed by Borue & Orszag (1996) is 
used as a specific example. 

The ensemble average is typically used for theoretical studies, the tem- 
poral average for experimental data, and the spatial average for data from DNS; 
thus all three are employed here. Ensemble, time and space averages are not 
interchangeable because the averages commute differently with differential op- 
erators within the dynamical equations. For the homogeneous case and infinite 
averaging volume, the spatially averaged equation ( |25| ) and the ensemble aver- 
aged equation ( [l7|) reduce to the same form, and similarly for the temporally 
averaged equation for the stationary case and infinite averaging time. 

Ongoing interest in turbulence intermittency includes accurate evalua- 
tion of inertial-range exponents of structure functions, for which purpose precise 
definition of an observed inertial range is needed. The third-order structure 
function can serve this purpose because it has a well-known inertial-range power 
law and the 4/5 coefficient (Kolmogorov 1941). Deviations from the 4/5 coeffi- 
cient are observed in experiments (Anselmet, Gagne & Hopfinger 1984, observe 
values greater than 4/5 in duct flow and less than 4/5 in jet flow despite the 
Reynolds numbers being nearly equal; Mydlarski & Warhaft 1996 and Lindborg 
1999 obtain values less than 4/5, etc.); this casts doubt on the precision with 
which measured exponents apply to the intermittency phenomenon (Sreenivasan 
& Dhruva, 1998). The equations derived here, when evaluated with data, can 
reveal all effects contributing to the deviation from Kolmogorov's 4/5 law and 
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his equation. 

The plan of the paper is to develop the mathematics in §2 and §3; §2 
contains necessary definitions and unaveraged equations; §3.1-3.3 contain the 
definition of averaging operations and their application to produce averaged 
equations. §3.4 contains the simplifications for the case of spatially periodic 
DNS. §3.5 defines the sphere average in r space and its associated orienta- 
tion average and relates these to the work of Kolmogorov (1962) and Obukhov 
(1962); Kolmogorov's equation is derived in §3.6 as a useful point of reference. 
Discussion of opportunities that these equations present for future investigations 
is in §4. 

2 Exact unaveraged two-point equations 

The equations given here relate two-point quantities and are obtained from 
the Navier-Stokes equations and incompressibility. The two spatial points are 
denoted x and x'; they are independent variables. They have no relative motion; 
e.g., anemometers at x and x' are fixed relative to one another. To be concise, 
velocities are denoted U{ = itj(x, t), u\ — u,(x',i), energy dissipation rates by 
e = e(x, t), e' = e(x', t) ,etc. p is the pressure divided by the density (density 
is constant), v is kinematic viscosity, and d denotes partial differentiation with 
respect to its subscript variable. Summation is implied by repeated Roman 
indices; e.g., d Xn d Xn is the Laplacian operator. For brevity, define: 

dij = (ui - u'j) (uj - u'j) , (1) 
d ijn = (ui - u'j) (uj - u'j) (u n - u' n ) , (2) 

T i3 = (d Xi p - d x > t p') (uj - u'j) + (d X] p - d^p'^j (Ui - u'i) , (3) 

eij = (d Xn Ui) (d Xn Uj) + (d x ' n u'i) (d x > n u'j) , (4) 

F l3n = {u t -u' i ){u J -u' j ) + ( 5 ) 

Change variables from x and x' to the sum and difference independent variables: 

X = (x + x') /2 and r = x — x', and define r = |r| . 
The derivatives dx t and d Ti are related to d Xi and d x > by 

d Xl = d ri + ^d Xz , d x >. = -d rt + ^d Xl , d Xi = d Xz +d x ^ , d ri = ^ (d Xi - d x >^) . 

(6) 

For any functions /(x, t) and g{x.' ,t), <J6|) gives 

d n [/(x, t) ± 5 (x', t)} = 3 Xi [/(x, t) T 5 (x', t)} j2. (7) 
Use of (^|) in (||) and in the trace of (Q) and rearranging terms gives 

Tij =-2(p- p) (sij - s'ij) + d Xi [(p - p) {uj - u'j)] + d Xj [{p - p) (ui - u'i)] , 

(8) 
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en = v- l {e + e')+d Xn d Xn {p+p'), (9) 

where Sjj = (d Xi uj + d Xj ui) /2 , and e = 2vsijSij. (10) 

To obtain (||) we used Poisson's equation, d x d Xn p = —d Xi Ujd Xj Ui. Incom- 
pressibility requires su = 0; thus, the trace of (B) is 

T ii = 2d Xi [(j)- P ')(u i -u'i)}. (11) 

Of course, all quantities above are local and instantaneous quantities. 



2.1 Use of the Navier— Stokes equation 

The Navier-Stokes equation and incompressibility give, 

d t dij+d Xn F ij n + d rn dijn = -Tij + 2v (d rn d rn d i: j + jd Xn d Xn d i: j ~ e ij^J ■ ( 12 ) 

As a check, one sees that ( fL2| ) can be obtained from equation (2.9) of Hill (2001). 
The trace of (|l^) and substitution of (||) and ( pd| ) gives 

dtdu + d Xn F ri n + d Tn d lin = 2vd rn d rn da -2(e + e') + w, (13) 

where 

w = -28 Xi [(p - p') { Ul - u'i)} + ^-d Xn d Xu d u - 2vd Xn d Xn (p + p') . (14) 



The first term in ( |14[) is — Tu fro m (|ll| ) and the last term in ( |14[) arises from 
en in (^|) ; the disparate terms in (|l4[) are given the symbol w for subsequent 
convenience and brevity. The limit r — > applied to (|13| ) recovers the definition 
of e in (fol). 



2.2 Exact second-order incompressibility relationships 

Because x and x' are independent variables, d Xi u'j = 0, and d x >_Uj = 0. 
Then, incompressibility gives: d Xn u n = 0, d Xn u' n = 0, d rn u n = 0, d rn u' n = 0, 
so <9x„ (u„ — u' n ) = 0, and d rn {u n — u' n ) = 0. The combined use of those 
incompressibility relations and (|7|) gives 

d r n [(uj - u'j) (u n - u' n )] = d Xn [(uj + Uj) (u n - u' n )] /2, (15) 
d rj d rn [(uj - vfj) (u n - u' n )] = d Xj d Xn [(uj + u'j) (u n + u' n )] /A. (16) 



3 Exact averaged two-point equations 

3.1 Ensemble average: exact equations 

The ensemble average is defined at each point (X, r, t) as the arithmetical 
average over the ensemble. We denote the ensemble average by angle brack- 
ets and subscript E: {o) E . Because the ensemble averaging operation is a 
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summation, it commutes with differential operators, the average of (|T^) is 
dt (dij) E + dx n (Fij tl ) E + 9 rn (dij n ) E 
= - ( T y ) E + 2v ( d r n d rn (dij ) E + 2 d x n d Xn (d^ ) E - (e i:i ) E ^j . (17) 



The argument list for each tensor in ( |17j ) is (X, r,i); the ensemble average does 
not eliminate dependence on any independent variable. The average of (|l3[) is 



d t {da) E + d Xn (F « n ) E + d Trl (d im ) E = 2vd Tn d rn (d u ) E - 2 (e + e') E + (w) E . 

(18) 

Exact incompressibility conditions on the second-order velocity struc- 
ture function are given by the ensemble averages of © and @: 

dr n (d 3 n) E - d Xn {( Uj + u' 3 ){u n -u' n )) E /2, (19) 
= dx 3 d Xn {{u J +u' j ){u n + u' n )) E /A. (20) 

3.2 Temporal average: exact equations 

Because nearly continuous temporal sampling is typical, we represent the 
temporal average by an integral, but all results are valid for the sum of discrete 
points as well. The temporal average is most useful when the turbulence is 
nearly statistically stationary. Let to be the start time of the temporal average 
of duration T. The operator effecting the temporal average of any quantity Q 
is denoted by (o) T , which has argument list (X,r,to,T); that is, 

(Q) T = ±[ t0+T Q(X,r,t)dt. (21) 
1 J to 



The argument list (X, rjo,T) is suppressed. The temporal average of (pL 2|— |l 6|) 
gives equations that are the same form as (p^7|-p0|) with one exception. The 



exception is that dt does not commute with the integral operator (21) such that 
(9tdij) T appears, whereas dt (dij) E appears in (|l7|), and similarly for the trace 
equation (|I"g|). Because data are taken at x and x' in the rest frame of the 
anemometers, and dt is the time derivative for that reference frame, it follows 
that 



(d t dij) T = 1 f ° d t dijdt = [dij (X, r, < + T) - dij (X, r, to)] /T. (22) 

1 Jto 



This shows that it is easy to evaluate (d t dij) T using experimental data because 
only the first (at t — to) and last (at t = to + T) data in the time series are 
used. One can make (dtdij) T as small as one desires by allowing T to be very 
large provided that d^ (X, r, to + T) does not differ greatly from dij (X, r, to). 
This is aided by judicious choice of to and to+T for the stationary case, but is 
not possible in all cases. 
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3.3 Spatial average: exact equations 

Because nearly continuous spatial sampling is typical of DNS, we represent 
the spatial average by an integral, but all results can be generalized to the case 
of a sum of discrete points. Let the spatial average be over a region R in X- 
space. The spatial average of any quantity Q is denoted by (Q) R which has 
argument list (r, f,R); that is, 

(Q)m^^J J jQ(K,r,t)dX, (23) 

where V is the volume of the space region R. The argument list (r, t, R) 
is suppressed. The spatial average commutes with r and t differential and 
integral operations, and with ensemble, time, and r-space averages, but not 
with dx n - Given any vector the divergence theorem relates the volume 
average of dx n qn to the surface average; that is, 

(d Xn q n ) R = ±- [ [ [ d Xn q n dX = ~ f ~ / / N n q n ds] = £ / q n , (24) 



V J J J " V \S J J / V _ 

where S is the surface area bounding R, dS is the differential of surface area, 
and N n is the unit vector oriented outward and normal to the surface. As 
seen on the right-hand side of (24), we adopt, for brevity, the integral-operator 
notation 

' 1 N n dS. 



/x n S 
The spatial average of ( |l2] ) is 

S 



d t (dij) R + — J Fijn + 9 rn (dijn) K 

= - ( t ^)k + 2v (^r n drn (dij) M + ~ d Xn dij - (ey) R ^ . (25) 
The spatial average of (|l^) is 
d t (dii) M + ^ <j> F un+dr n (dun) R = 2vdr n d Vn (di<) R -2 (e + e') R + (w) 
S 



(26) 



where (w) K _ 

v Jx„ 



-2 (p - p') (u n - u' n ) + -dx n dij - 2ud Xn (p + p) 



The spatial average of the incompressibility condition ( |15|) is 

d r„ (djn) R = M n) («j + u j) » ( 27 ) 

which is, on the right-hand side, a surface flux of a quantity that depends on 
large-scale structures in the flow. Similarly, (Qq) gives 



d rj d rn (d jn ) u = d Xj [K + <) («j + «j-)] ■ (28) 



6 



3.4 Spatial average: DNS with periodic boundary condi- 
tions 



The spatial average is particularly relevant to DNS. DNS that is used to 
investigate turbulence at small scales often has periodic boundary conditions. 
For such DNS, consider the spatial average over the entire DNS domain. Con- 
tributions to ^ x q n from opposite sides of the averaging volume cancel for that 
case such that ^ x q n = and therefore (dx n Qn) R = 0. In (|25| ) we then have 
ix Fijn — and jf x dx n dij = 0. In j26| ) we have ^ x Fu n = and {w) R = 0. 
In (p7|), the right-hand side vanishes so that d Tn (dj n ) R = 0. Thus, in the DNS 
case described above, we have 

d t (c%) K + d Tn (d ijn ) R = - (Tij) M + 2v (d rn d rn (dij) M - (ey) M ) , (29) 
d t {du) R + d rn (d iin ) u = 2vd rn d rn {du) R - 2 (e + e') R , (30) 

dr n (djn) R = 0, and 5 r „ (e,-„) R = 0. (31) 
Proof of d Tn (e 3n ) R = is given in Hill (2002). 



Performing the r-space divergence of fl29|) and using (31), we have 



drjdr n (dij n ) R — d rj (Tij) R . (32) 

This exact result is analogous to the asymptotic result in Frisch (1995), Lindborg 
(1996), and Hill (1997). 

Using the Taylor series of e and e' around the point X, Hill (2002) 
obtains the following exact result for the periodic DNS case considered 

-2(6 + 6% = -4(e(X,i)) R , and - (e^ = -4v([(d Xn Ui) (9 In «j)] x=x L , 

(33) 

where the subscript x = X means that the derivatives are evaluated at the point 
X. An important aspect of ( |33| ) is that the right-hand sides depend only on t. 
Of course, none of (p9|-|32"|) depends on X because of the spatial average over X. 

No approximations have been used to obtain the above equations for 
the spatially periodic DNS case considered. 

3.5 Averages over an r-space sphere 

The energy dissipation rate averaged over a sphere in r-space has been 
a recurrent theme in small-scale similarity theories since its introduction by 
Obukhov (1962) and Kolmogorov (1962). By averaging our equations for the 
trace, we can, for the first time, produce an exact dynamical equation containing 
the sphere-averaged energy dissipation rate. The volume average over an r- 
space sphere of radius rg of a quantity Q is denoted by 



(QWere = (^|/3) 'J J Jq (X, r, t) dr. (34) 



|r| < rs 
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The orientation average over the surface of the r-space sphere of radius r$ of a 
vector q n (X, r, t) is denoted by the following integral-operator notation: 

j> q n = (AnrQ- 1 J J ?fq n (X, r, t) ds, (35) 

|r| = r s 

where ds is the differential of surface area, and r n /r is the unit vector ori- 
ented outward and normal to the surface of the r-space sphere. Note that 
(47rr|) ds = dfi / 47r where dtl is the differential of solid angle from the sphere's 
center. Both {Q) r _ sp i lcre an d § r q n have the argument list (X,rs,i), which is 
suppressed. The divergence theorem is 

(dr„q n ) r . Bphcic = (3/r s )^ qn- (36) 

Because r, X, and t are independent variables, the r-space volume and ori- 
entation averages commute with time and X-space averages and with X- and 
^differential operators, and, of course, with the ensemble average as well. For 
instance, 

{d t (dii)m)r-sphere = dt (^"Mr-sphorc = ( (^W™) r-sphere) = 9 t (( rf n) 



r-sphere / ' e ^ C ' 



The r-sphere average (g4j) can operate on all of the above structure- 
function equations; it can operate on unaveraged equations ( |l2|) and (|l^) as well. 
These equations have terms of the form d Tn q n \ e.g., q n = (c% n ) K , d Tn (du) Kl 
(dun) ei (dijn) T , d rn (du) T , etc. By means of ([36]), the volume average in r-space 
of any term of the form d Tn q n produces the orientation average of q n within the 
subject equation. After operating on ( |l8| ) with the volume average in r-space 

(H]), the term — 2(e + e') E in that equation produces — 2 ^(e + e') r _ sphcrc ) . 

Now, ^(^ + s') r-sphere) ^ s tne sphere-averaged energy dissipation rate de- 
fined in the third equations of both Obukhov (1962) and Kolmogorov (1962). 

The result of the r-space sphere average of any of our equations will be 
clear from operating on (^). The average of ( J30|) over a sphere in r-space of 
radius rs and multiplication by rs/3 and use of (|33|) gives 



r f9t((du) r _ sphc , c ) R +f (d iin ) K = 2vf d rn (du} M -^((e) r _ spheIe ) R . (37) 



The terms have argument list (rs,t), but ( (£} r _ sp h C rc) depends only on t. Of 

course, none of the quantities in (|37|) depends on X because of the X-space 
average. Despite its simplicity, ( |37| ) has been obtained without approximations 
for the freely-decaying spatially-periodic DNS case considered; ( |37| ) applies to 
inhomogeneous and anisotropic DNS that have periodic boundary conditions. 
Nie & Tanveer (1999) define a structure function S3 using time, space, and 
solid-angle averages acting on du n , and consider the asymptotic inertial range 
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case to obtain that £3 = — (4/3) er without use of local isotropy. An analogous 
result can be obtained by applying inertial-range asymptotics to (|37|); namely, 
neglect the time-derivative term on the basis of local stationarity and neglect 
the term proportional to v. 



3.6 Kolmogorov's equation derived from (|37|) 

Most readers are familiar with Kolmogorov's (1941) famous equation that 
is valid for locally isotropic turbulence. A useful point of reference is to derive 



it from (37). This helps elucidate (J37|). An index 1 denotes projection in the 
direction of r, and indices 2 and 3 denote orthogonal directions perpendicular 
to r. For locally isotropic turbulence we recall that the only nonzero compo- 
nents of (d ijn ) R are (d m ) K , (d 22 i) R = (^33i} R , and of (dy) B are (du) R , and 
(<^22}r = (<^33)r- These components depend only on r such that there is no 
distinction in an r-space sphere average between rs and r; thus, we simplify 
the notation by replacing rs with r. The isotopic-tensor formula for (dij n ) R 
gives (diin) m = (r n /r) ((dm) R + 2 (cfei)^) = (r n /r) (dui) m substitution of 
which into @ gives ^ (d iin ) M = (r„/r) {d iin ) R = (r n /r)(r n /r) {d ul ) R = 
{dui) m - Since (9 r „r) = (r„/r), we have § tn <9 r „ {d u ) R = (r n /r) (d rn r) d r («fe) R = 
d r (du) R . Kolmogorov (1941) considered the locally stationary case such that 
he neglected the time-derivative term, thus we also neglect that term to obtain 
from <M) that 

4 

(din) M = 2vd r {du) u - - (e) K r. (38) 



Alternatively, we can time average (pTj); then the time derivative can be ne- 
glected with the weaker condition noted with respect to the smallness of (p2|); 



then ((dm) R ) T = 2vd r ((<fa) R ) T — | ((e) R ) T r. For simplicity of notation, wc 
continue with (|3|). To eliminate (e£22) R and ((d22i)R) from the expressions 
(du) M = (dn) R + 2(d 22 ) R and (d U i) M = (dm) H + 2 (d 22 i) R , we use the in- 
compressibility conditions |9 r (dn} R + (dn) R — (d22) R — 0, and rd r (dm) R + 
(^iii)k — 6 (^22i)r = 0, which are valid for local isotropy (Hill 1997), and were 
also used by Kolmogorov (1941). Then ( p8| ) becomes, after multiplying by 3r _1 , 
d r (dm) M + 7 ( d m) R = § v [ d r (dn) R + 7~(dii) R ] - 4 (e) R ; this is Kolmogorov's 
(1941) third equation. After multiplication by r 4 and integrating from to r, 
we have Kolmogorov's (1941) equation 

4 

(dm) R = Qvd r {dn) R - - (e) R r. (39) 

Kolmogorov's inertial-range 4/5 law and the viscous-range law follow immedi- 
ately from @. 
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4 Examples of opportunities for using the exact 
equations 



4.1 Mitigating anisotropy to check asymptotic laws 

Consider homogeneous, anisotropic turbulence. Homogeneity causes dx n 
operating on a statistic to vanish (Hill 2001), so dx n (F ij n ) E and dx n dx n (dij) E 
vanish from (|l~7|), but (rij) E becomes — 2((p — p') (s^ — s'ij)) E (see ||), which 
does not vanish. Under the more restrictive assumption of local isotropy, 
( T y) e = ^ (-^ m 1997) such that the entire nonzero value of (Tij) E is a source 
of anisotropy in (|T7|). For the locally stationary case, the anisotropy quantified 
by (Tjj) E is approximately balanced by that from the term d Trl (dij n ) E in (O) 
(Hill 1997, and exactly so for the stationary case). In contrast consider (|18|). 
Homogeneity causes dx n F u n and w to vanish from (p"8|); equivalently, (tu) e 
is absent from ( |l8| ) because incompressibility gives sa — s' u — 0. Therefore, 
for the homogeneous, anisotropic case, an important source of anisotropy of 
d rn (dij n ) E , namely (ry) E , is absent from d Tn (da n ) E . It therefore seems that 
(dm) E will more accurately show the asymptotic inertial-range power law than 
does (dm) E (or (d 2 2i) E or (ds3i) E ). This result for the homogeneous case 
extends to the locally homogeneous case as follows: For inhomogeneous turbu- 
lence, the nonvanishing part of (to) E , i.e., (tu) e = 2dx t ((p — p') (ui — u' i )) E 
[see (|ll|)1 is expected to approach zero rapidly as r decreases for two reasons. 
First, ((p — p') (ui — u' i )) E vanishes on the basis of local isotropy. Second, 
the operator dxi causes &Xi ((p —p') — u 'i)) e to van i sn on the basis of lo- 
cal homogeneity. From (|l4|), (w) E contains the terms vdx n dx n (du) E /2 and 
—2vdx n dx n + because of the operator dx n dx n , these terms vanish on 
the basis of local homogeneity. Thus, all terms in (w) E are negligible for locally 
homogeneous turbulence. By performing the trace, it appears that anisotropy 
has been significantly reduced in ( |l8| ) relative to in ( p~7| ) for the high-Reynolds- 
number, locally homogeneous case such that the above hypothesis is extended 
to locally homogeneous turbulence. The hypothesis should be checked by com- 
parison with anisotropic DNS. Evaluation of all terms in ( |l7| ) and (|l^) are the 
basis for such an investigation. The above discussion holds for temporal and 
spatial averages as well. 

To determine scaling properties of the third-order structure function, 
past theory has used the isotropic formulas. One can use an equation like ( p"8| ) 
or its temporal-average analogue without an assumption about the symmetry 
properties (e.g., isotropic) of the structure functions by means of the sphere av- 
erage in r-space. Without approximation, the r-space sphere average produces 
the orientation-averaged third-order structure function. It would seem that the 
orientation average mitigates anisotropy effects. Thus, the orientation average 
§ Tn (divn) E ( or § rn (d iin ) T , or § r (d im ) R ) is expected to best exhibit scaling 
properties of locally isotropic turbulence, such as the inertial-range power law 
with the 4/3 coefficient in (|37|). This hypothesis should be checked by compar- 
ison with anisotropic DNS. 



10 



Consider the stationary, homogeneous case. From (|18|), we are certain 
that d rn (dun) E — 2vd rn d rn (du) E scales with (e + e') E because 

(dm ( d an) E ~ 2 ^ r „d r?l (d u ) E ) / (e + e') E = -2; 

thereby insuring K41 scaling of d rn {du n ) E — 1vd Tn d Tn (dy) E despite anisotropy. 
In contrast, (|T7]) insures that scaling only if local isotropy is invoked. Anisotropic 
DNS can be used to check whether or not K41 scaling is improved by performing 
the trace. 



4.2 Tests using DNS and experimental data 

The spatially periodic DNS case leads to especially simple equations. It 
seems that (p9[-|30|) offer an ideal opportunity to evaluate the contribution of 
dt (dij) M for freely decaying turbulence, and of (Ty ) R for anisotropic turbulence, 



as well as the balance of the off-diagonal components of (29). Because we have 
not introduced a force generating the turbulence and because every point in 
the flow enters into the X-space average, the DNS must be freely decaying. 
As shown in Hill (2002), it is straightforward to include forces in the equations. 
New experimental methods (Su & Dahm 1996) and DNS can completely evaluate 
terms in the exact structure-function equations. 



4.3 Effect of inhomogeneity on incompressibility condi- 
tions 

Exact incompressibility relationships (|l9|-|20|) are obtained that can be used 
to quantify the nonzero value of d Tn (dj n ) E (or of d Tn (dj n ) T , or of d Tn (dj n ) R ) 
caused by inhomogeneity. If inhomogeneity is only in the streamwise (say 1- 
axis) direction, then the time average gives d Tn (dj n ) T — dx x ((uj + u'j) (ui — u i)) 
which can be evaluated using anemometers. As r — > 0, (Um becomes the sec- 
ond derivative with respect to measurement location of the velocity variance 
and therefore clearly depends on flow inhomogeneity. 

4.4 Quantifying effects of inhomogeneity and anisotropy 
on scaling exponents 

Sreenivasan & Dhruva (1998) note that one could determine scaling ex- 
ponents with greater confidence if one has a theory that exhibits not only the 
asymptotic power law but also the trend toward the power law. Such a theory 
must require difficult measurements or DNS to evaluate such trends. The equa- 
tions given here are the required theory for the third-order structure function, 
given that data must be used to evaluate the equations in a manner analogous 
to previously cited evaluations. In fact, it is not possible that exact equations 
do not contain the physical effects discussed by Sreenivasan & Dhruva (1998). 
They discuss the fact that there is correlation of velocity increments with large- 
scale velocity in inhomogeneous turbulence, even for very large Reynolds num- 
bers and r in the inertial range, but not so in isotropic turbulence. Our term 
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dx„ (F un) E — dx n — u '| 2 (u n + u n) /2 ) in (^8|), and all such analogous 
terms in the other equations, explicitly contains such correlation, and the bal- 
ance of the equations imparts that correlation effect to the other statistics; all 
such terms do vanish for isotropic turbulence. They also discuss the usefulness 
of graphing all 3 terms in (39) to discern the onset of the dissipation-range. 
Our equations are exact there too. 



4.5 Quantifying effects of large-scale structure on small- 
scale structure 

Experimenters remove the mean from an anemometer's signal before cal- 
culating structure functions from the velocity fluctuations, whereas the exact 
dynamical equations contain statistics of the full velocity field. Hill (2002) 
applied the Reynolds decomposition to the above exact dynamical equations, 
and used inertial-range and viscous-range asymptotics to determine the approx- 
imate dynamical equations pertaining to statistics of fluctuations as well as all 
approximations that are required to obtain the approximate equations. The 
Reynolds decomposition produces terms that quantify the effect of the large- 
scale structure of turbulence on the small scales. For example, dx n (F ijn) E 
produces a generalization of the advective term discovered by Lindborg (1999). 
Hill (2002) contrasts the various definitions of local homogeneity and points out 
that the only definition that simplifies dynamical equations is that from Hill 
(2001). 
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